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^T \ Abstract 

o 

^^ . Supersymmetric extensions of the Standard Model with small R-parity and lepton 

. ■ number violating couplings are naturally consistent with primordial nucleosynthe- 

r—-j ■ sis, thermal leptogenesis and gravitino dark matter. We consider supergravity mod- 

els with universal boundary conditions at the grand unification scale, and scalar 
r- lepton or bino-like neutralino as next-to-lightest superparticle (NLSP). Recent 
Fermi-LAT data on the isotropic diffuse gamma-ray flux yield a lower bound on 
the gravitino lifetime. Comparing two-body gravitino and neutralino decays we 



find a lower bound on a neutralino NLSP decay length, cr v o > 30 cm. Together 
with gravitino and neutralino masses one obtains a microscopic determination of 
the Planck mass. For a t-NLSP there exists no model-independent lower bound 
on the decay length. Here the strongest bound comes from the requirement that 
the cosmological baryon asymmetry is not washed out, which yields ct~ > 4 mm. 
However, without fine-tuning of parameters, one finds much larger decay lengths. 
For typical masses, m 3 y 2 ~ 100 GeV and mNLSP ~ 150 GeV, the discovery of a 
photon line with an intensity close to the Fermi-LAT limit would imply a decay 
length ctnlsp of several hundred meters, which can be measured at the LHC. 



1 Introduction 

Locally supersymmetric extensions of the Standard Model predict the existence of the 
gravitino, the gauge fermion of supergravity pQ. For some patterns of supersymmetry 
breaking, the gravitino is the lightest superparticle (LSP), and therefore a natural dark 
matter candidate [2J. Heavy unstable gravitinos may cause the 'gravitino problem' [3L 
[5] for large reheating temperatures in the early universe. This is the case for thermal 
leptogenesis [6] , where gravitino dark matter has become an attractive alternative [7] to 
the standard WIMP scenario [5]. 

Recently, it has been shown that models with small R-parity and lepton number 
breaking naturally yield a consistent cosmology incorporating primordial nucleosynthesis, 
leptogenesis and gravitino dark matter [9J. The gravitino is no longer stable, but its 
decays into Standard Model (SM) particles are doubly suppressed by the Planck mass 
and the small R-parity breaking parameter. Hence, its lifetime can exceed the age of the 
Universe by many orders of magnitude, and the gravitino remains a viable dark matter 
candidate [TU] . 

Gravitino decays lead to characteristic signatures in high-energy cosmic rays, in par- 
ticular to a diffuse gamma-ray flux [9TTTT]. The recent search of the Fermi-LAT collabora- 
tion for monochromatic photon lines [TS] and the measurement of the diffuse gamma-ray 
flux up to photon energies of 100 GeV [19] severely constrain possible signals from de- 
caying dark matter. In this paper we study the implications of this data for the decays 
of the next-to-lightest superparticle (NLSP) at the LHC, extending the estimates in [9J. 

We shall restrict our analysis to the simplest class of supergravity models with univer- 
sal boundary conditions at the Grand Unification (GUT) scale, which lead to neutralino 
or f-NLSP. Electroweak precision tests, thermal leptogenesis and gravitino dark mat- 
ter together allow gravitino and NLSP masses in the range m^/2 = 10 . . . 500 GeV and 
ttjnlsp = 100 . . . 500 GeV [20] ■ Following [9], the breaking of R-parity is tied to the break- 
ing of lepton number, which leads to a model with bilinear R-parity breaking [2T|I22] . The 
soft supersymmetry breaking terms are characteristic for gravity or gaugino mediation. 

In order to establish the connection between the gamma-ray flux from gravitino de- 
cays and NLSP decays, one needs R-parity breaking matrix elements of neutral, charged 
and supercurrents. For the considered supergravity models we are able to obtain these 
matrix elements to good approximation analytically. This makes our results for the NLSP 
decay lengths rather transparent. As we shall see, the lower bound on the neutralino de- 
cay length is a direct consequence of the Fermi-LAT constraints on decaying dark matter. 
On the other hand, the lower bound on the f-decay length is determined by the cosmo- 
logical bounds on R-parity breaking couplings, which follow from the requirement that 



the baryon asymmetry is not washed out 

This paper is organized as follows. In Section 2 we discuss the general Lagrangian 
for R-parity breaking in a basis of scalar SU(2) doublets where all bilinear mixing 
terms vanish. This leads to new Yukawa and gaugino couplings, some of which are 
proportional to the up-quark Yukawa couplings. Section 3 deals with the various 
supersymmetry, R-parity and lepton number breaking terms in the Lagrangian and the 
relations among them due to a U(l) flavour symmetry of the considered model. The 
needed R-parity breaking matrix elements of neutral, charged and supercurrent are 
analytically calculated in Section 4, based on the diagonalization of the mass matrices 
which is discussed in detail in the appendix. The main results of the paper, the bounds 
on the NLSP decay lengths and the partial decay widths, are described in Section 5, 
followed by our conclusions in Section 6. 



2 Bilinear R-Parity Breaking 



Supersymmetric extensions of the Standard Model with bilinear R-parity breaking con- 
tain mass mixing terms between lepton and Higgs fields in the superpotential 
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AW = fuHJi (1) 

as well as the scalar potential induced by supersymmetry breaking, 

- AC = BiHJi + m 2 J\H d + h.c. . (2) 

These mixing terms, together with the R-parity conserving superpotential 

W = fiH u H d + h u l3 q t u c 3 H u + h%d c % q 3 E d + h%ke)H d , (3) 

the scalar mass terms 

-C M =m 2 u HlH u + m 2 d H ] d H d + (BH u H d + h.c.) 

+ rnffiiU + ml^i^i + ™ 2 qi qUi + m 2 ui ufu\ + m\dfd\ , (4) 

and the standard SU(3) x SU{2) x U(l)y gauge interactions define the supersymmetric 
standard model with bilinear R-parity breaking. Note that the Higgs mass terms rr? u and 
m 2 , contain the contributions from the superpotential ([3]) and the soft supersymmetry 
breaking terms. For simplicity, we have assumed flavour diagonal mass matrices in (J4J). 



^^Our notation for Higgs and matter superfields, scalars and left-handed fermions reads: H u 
(H u ,h u ), h = (k,U) etc. 



For a generic choice of parameters the electroweak symmetry is broken by vacuum 
expectation values (VEVs) of all scalar 577(2) doublets, 

{H° u )=v u , (H° d )=v d , <*,> = «,, (5) 



witltl 

v u _ . n ~ _ Vi _ Bj tan P-mf d -fifi* 

— = Tan p , 6i = — — i — 2 7T3 — > \°) 

Vd v d mfi + 5 m| cos 2p 

where higher order terms in the R-parity breaking parameters have been neglected. 

It is convenient to discuss the predictions of the model in a basis of SU(2) doublets 
where the mass mixings /ij, 5, and m 2 d in Eqs. ([T]) and ([2]) are traded for R-parity 
breaking Yukawa coulings. This can easily be achieved by field redefinitions. First one 
rotates the superfields H d and /,, 

Hd = H d — e i l i , k = l { + 6iH d , €i — — . (7) 

Then the bilinear term ([I]) vanishes for the new fields, i.e., /^ = 0, and one obtains 
instead the cubic R-parity violating terms 

AW = \\^ % e% + \[ jk dtvA , (8) 

where 

Kjk = -hij€ k + h e kj ei , X ijk = -h tj e k . (9) 

The new R-parity breaking mass mixings are given by 

B' i = B i -Be i , mf d = m 2 id + ei{ml - m 2 d ) . (10) 

The corrections for R-parity conserving mass terms are negligable. 

In a second step one can perform a non-supersymmetric rotation among all scalar 
SU(2) doublets, 

H d = H^-e^, eH* u = eH>: - tf! , ? = If + ^H" d + &HZ , (11) 

where e is the usual S77(2) matrix, e = ir 2 . Choosing 

B'.B + mlJml-ml) 
l ~ {ml-jnl){ml-m 2 d )-B 2 > l ) 

„ = B\ (m 2 u - m 2 d ) + Bmf d 
1 (™l - m 2 ) (m 2 H - ml) - B 2 ' l ) 



2 Note that our result for Zi = Vi/vd holds at all renormalization scales, contrary to different expres- 
sions used in the literature. 



the H u li and VHd mixing terms vanish in the new basis of doublets. According to (JH]) 
also the scalar lepton VEVs (pi) vanish in this basis. 

It is straightforward to work out the R-parity violating Yukawa couplings which are 
induced by the rotation (ITTj) . We are particularly interested in the terms containing one 
light superparticle, i.e, a scalar lepton, bino or wino. The corresponding couplings read, 
after dropping prime and double-prime superscripts on all fieldgj, 
1 
2' 



■A£ d -Xijkli&jlk + ^ijkd'iQjh + \jkhejlk + ^ijk ( li u 'j £ lk 



+ htMH d + e'leH:)e c j k 



<i 



,/ 



- JL&Hl - 4Hle)kb + -j=(4Hl - ^HleVW + h.c. , (14) 

where the Yukawa couplings are given by 

Aijfc = ~h%e k + h e kj e { , \' ijk = -h%{e k + e' k ) , (15) 

Kjk = -hij(e k + e k ) + h e kj €i , X ijk = h^e k . (16) 

Since the field transformations are non-supersymmetric, the couplings \j k and \j k are 
no longer equal as in Eq. §§§. Furthermore, a new coupling of right-handed up-quarks, 
\'ij k1 has been generated. 

After electroweak symmetry breaking one obtains new mass mixings between higgsi- 
nos, gauginos and leptons, 

-A£ M D m^— e^hd-mzSwQvib + m z Cw£viW 3 + h.c. , (17) 

J cp J 

where we have defined 

0= ffoi±jk v = ^Tv>, i = ^ s S, (18) 

V V V d Cfi 

y/g 2 + g' 2 v gf_ 

V2 ' Sw ~ Vg^+7 1 



m ij = he ijVd , m z = j= s w = ,_,, : r> = -y/l - 4 . (19) 



Given the Yukawa couplings hfj, hfj and /i?-, the Lagrangian (TT41) predicts 108 Im- 
parity breaking Yukawa couplings in terms of 9 independent parameters which may be 
chosen as 

IH , Bi , m 2 id or e 4 , e^ , e'[ . (20) 

These parameters determine lepton-gaugino mass mixings, lepton-slepton and quark- 
slepton Yukawa couplings, and therefore the low-energy phenomenology. The values of 
these parameters depend on the pattern of supersymmetry breaking and the flavour 
structure of the supersymmetric standard model. 



3 Our notation for gauge fields, field strengths and left-handed gauginos reads: i? M , B^, b etc. 
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3 Spontaneous R-parity breaking 

Let us now compute the parameters e», e' t and e'{ in a specific example where the sponta- 
neous breaking of R-parity is related to the spontaneous breaking of B-L, the difference 
of baryon and lepton number [9] . 

We consider a supersymmetric extension of the standard model with symmetry group 

G = SU(3) x SU(2) x U(1) Y x U(1) b ~l x U(1) r . (21) 

In addition to three quark lepton generations and the Higgs fields H u and H d the model 
contains three right-handed neutrinos uf, two non-Abelian singlets N c and N, which 
transform as v c and its complex conjugate, respectively, and three gauge singlets X, $ 
and Z. The part of the superpotential responsible for neutrino masses has the usual form 

w v = KM H - + -^ h l v >l N2 > ( 22 ) 

where Mp = 2.4 x 10 18 GeV is the Planck mass. The expectation value of H u generates 
Dirac neutrino masses, whereas the expectation value of the singlet Higgs field N gen- 
erates the Majorana mass matrix of the right-handed neutrinos vf. The superpotential 
responsible for B-L breaking is chosen as 

W B -l = X(NN C - $ 2 ) , (23) 

where unknown Yukawa couplings have been set equal to one. $ plays the role of a 
spectator field, which will finally be replaced by its expectation value, (<&) = vb-l- 
Similarly, Z is a spectator field which breaks supersymmetry and U(1)r, (Z) = Fz09. 
The superpotential in Eqs. (|22|) and (|2"3"|) is the most general one consistent with the 
R-charges listed in Table [H up to nonrenormalizable terms which are irrelevant for our 
discussion. 

The expectation value of $ leads to the breaking of B — L, 

(N) = (N c ) = ($) = vb-l , (24) 

where the first equality is a consequence of the U(1)b-l D-term. This generates a 
Majorana mass matrix M for the right-handed neutrinos with three large eigenvalues 
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H u 


H d 


N 


N c 


$ 


X 


Z 


R 


1 











-2 


-1 
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Table 1: R-charges of matter fields $ — q, u c , e c , d c , I, v c , Higgs fields and gauge singlets. 



M 3 > M 2 > M\. If the largest eigenvalue of h n is 0(1), one has M 3 ~ v%_ L /M P . Inte- 
grating out the heavy Majorana neutrinos one obtains the familiar dimension-5 seesaw 
operator which yields the light neutrino masses. 

Since the field $ carries R-charge — 1, the VEV (<&) breaks R-parity, which is con- 
served by the VEV (Z) . Thus, the breaking of B — L is tied to the breaking of R-parity, 
which is then transmitted to the low-energy degrees of freedom via higher- dimensional 
operators in the superpotential and the Kahler potential. Bilinear R-parity breaking, as 
discussed in the previous section, is obtained from a correction to the Kahler potential, 

AK = -L (a i Z*&N c HJ i + a^Z^N^HJi 

IVi-p ^ 



- x htfZ&N c H u k + b'^ZQrfHuk 



M P 

+ Cl Z ] Z& N c l\H d + c'^Z^NHlHa) + h.c. . (25) 



Replacing the spectator fields Z and $, as well as N c and N by their expectation values, 
one obtains the correction to the superpotential 

AW = fnHJi , 

with 

^ = ^3(0^ + 0^3/26, e = ^=£~-^, (26) 

where m 3 / 2 = Fz/(-\/3Mp) is the gravitino mass. Note that G can be increased or de- 
creased by including appropriate Yukawa couplings in Eqs. (I2"2"j) and (1251) . The corre- 
sponding corrections to the scalar potential are given by 

- AC = BtHji + m 2 J\H d + h.c. , 

where 

Bi = 3(6; + 6-) m 3/2© , m 1d = 3(cj + Ci)ml /2 Q . (27) 

The corresponding R-parity conserving terms are generated by [2H] 

K D ^-&H u H d + J^&ZH u H d + h.c. , (28) 

IVlp IVlp 

which yields 

W D liH u H d , [i = y/3a m 3 / 2 , (29) 

-C D BH u H d + h.c. , B = 3b m 2 3/2 . (30) 

Higher dimensional operators yield further R-parity violating couplings between 
scalars and fermions. However, the cubic couplings allowed by the symmetries of our 



model are suppressed by one power of M-p compared to ordinary Yukawa couplings and 
cubic soft supersymmetry breaking terms. Note that the coefficients of the nonrenormal- 
izable operators are free parameters, which are only fixed in specific models of supersym- 
metry breaking. In particular, one may have /i 2 ,m 2 > m\/ 2 and hence a gravitino LSP. 
All parameters are defined at the GUT scale and have to be evolved to the electroweak 
scale by the renormalization group equations. 

The phenomenological viability of the model depends on the size of R-parity break- 
ing mass mixings and therefore on the scale vb-l of R-parity breaking as well as the 
parameters a\ . . .d i in Eq. (125J1 . Any model of flavour physics, which predicts Yukawa 
couplings, will generically also predict the parameters Oj . . . d { . As a typical example, we 
use a model [27] for quark and lepton mass hierarchies based on a Froggatt-Nielsen £7(1) 
flavour symmetry, which is consistent with thermal leptogenesis and all contraints from 
flavour changing processes |28j. 

The mass hierarchy is generated by the expectation value of a singlet field <fi with 
charge Q<f> = — 1 via nonrenormalizable interactions with a scale A = (<f))/i] > Aqut, 
t] ~ 0.06. The 77-dependence of Yukawa couplings and bilinear mixing terms for multiplets 
rpi with charges Qi is given by 

hij oc r] Ql+Q] , fa oc rp\ /3; oc 77 Q \ m 2 d ocr] Ql . (31) 

The charges Qi for quarks, leptons, Higgs fields and singlets are listed in Table [2J The 
neutrino mass scale m u ~ 0.01 eV implies for the heaviest right-handed neutrinos M 2 ~ 
M3 ~ 10 12 GeV. The corresponding scales for B — L breaking and R-parity breaking are 



v 2 
10 15 GeV , 6 = -^ ~ 10~ 6 . (32) 



M P 

For the small R-parity breaking considered in this paper the neutrino masses are domi- 
nated by the conventional seesaw contribution [PJ. 

The R-parity breaking parameters \x%i Bi and m 2 d strongly depend on the mechanism 
of supersymmetry breaking. In the example considered in this section all mass parameters 
are 0(7723/2), which corresponds to gravity or gaugino mediation. From Eqs. ( I2"6"l) .( j2"7j) 
and ( 13TT) one reads off 

^ = ari Q *m 3/2 , B t = bri Q >m 2 3/2 , m 2 id = cr] Q >m 2 3/2 Q , (33) 



tpi 


10 3 


10 2 


10i 


5s 


°2 


5^ 
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v\ 


H u 


H d 


$ 


X 


Z 


Qi 
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2 


1 


1 


2 








1 


















Table 2: Chiral U(l) charges. 10 4 = (q t , u 4 c , ef), 5 = {df, k), i = 1 . . . 3. 



with d,b,c = 0(1). Correspondingly, one obtains for e-parameters (cf. (IT21) . (IT3|) ) 

e . = a r) Qi Q , e[ = bri Q *Q , e'- d = cq Qi Q , (34) 

with a,b,c = 0(1). Our phenomenological analysis in Section 5.2 will be based on this 
parametrization of bilinear R-parity breaking. 

Depending on the mechanism of supersymmetry breaking, the R-parity breaking soft 
terms may vanish at the GUT scale [22], 

J B i (A GUT )=m2 d (A GUT ) = 0. (35) 

Non-zero values of these parameters at the electroweak scale are then induced by radiative 
corrections. The renormalization group equations for the bilinear R-parity breaking mass 
terms read (cf. [22], t = In A): 

16tt 2 ^ = 3/i 4 (h] k hf k - \g\ - gl) + n^hft - p (X ljk K 3 + 3X' kji h%) , (36) 



IGtt 2 ^ 
dt 


= 3B t (h%h% - \g\ - gt) + 6* ^M 1 + g\M^ 




+ B k h^h k * - B [Xijkhl* + 3X' kji h k *j , 


dt 


= *kji"'kj m d ~ rn jd' l jk' l ik ~ ^kji' l kj' m d + * l jk' l jk rn id 




+ 3h k ^h k ^m id + tn u X nki h nk — 3m u X nki h nk 




+ 2X k j i m lk X k j + 2X k ^ i h k ^m e j — vX k j i h k j , m dk — o\ k j i m q jh k j . 



(37) 



(38) 

In bilinear R-parity breaking, the R-parity violating Yukawa couplings vanish at the GUT 
scale. One-loop radiative corrections then yield for the soft terms at the electroweak scale 
(cf. Eqs. fl35D,([37D; €i = ^/fi): 

V* ( 6 „>2*S , «„2 /1/f \i A GUT 



^(A EW ) = -f^ (-9 ,2 M 1 + Qg 2 M 2 hi -^ , m 2 d (A EW ) = . (39) 

lb7i z \5 / A EW 

This illustrates that the bilinear R-parity breaking terms /z 2 , Bi and rn\ d are not neces- 
sarily of the same order of magnitude at the electroweak scale. 



4 Neutral, charged and supercurrents 

In Section 2 we have discussed the R-parity breaking Yukawa couplings in our model. For 
a phenomenological analysis we also need the couplings of the gauge fields, i.e., photon, 
W-bosons and gravitino, to charged and neutral matter, 

M c w M y/2 * a/2 m 2M P ^ v ; 



The corresponding currents read 



J, 



efi 



=w + 7 M w + - w 7 M w - ea^et + e^e* - h d ^h d + h u 7^u 



JztM =w + 7 M w + - w 'j^w + -va^Vi 



~.°ilpfii 






~Kivh d 



1 



1-0 



+ nKl»K ~ nKlvK 



s 2 J 



J? =v / 2 (vf-f^w - w + 7^w 3 ) + ucifjfii + h ( ff^h d + /^7/X , 

s» =\ [Y, Y] Y (bB up + w J wl p ) + .... 

The gravitino and the gauginos are now Majorana fermions, 

b = b + b c , w 1 = w 1 + w cI , 



(41) 

(42) 
(43) 
(44) 

(45) 



where the superscript c denotes charge conjugation. In Eqs. (l4"Tj) - (|4"4"1) we have only listed 
contributions to the currents which will be relevant in our phenomenological analysis. 

The R-parity breaking described in the previous section leads to mass mixings be- 
tween the neutralinos b, w 3 , /i°, h d with the neutrinos z/j, and the charginos w + , h+, w~, 
h d with the charged leptons e^, e^, respectively. The 7x7 neutralino mass matrix reads 
in the gauge eigenbasis 



M 



N 



( Mi 





m z sps w 


-m z c f5 s w 


-Q m zs w i 





M 2 


-m z s/3c w 


m z cpc w 


dm z c w 


m z S/3S w 


-m z spc w 





-P 





-m z cps w 


m z cpc w 


-\x 








\ -Cifn z s w 


d m zc w 








o ) 



(46) 



where we have neglected neutrino masses. Correspondingly, the 5x5 chargino mass 
matrix which connects the states (w~ , h d , e,) and (w + , h£, ef) is given by 



M u = 



I M 2 m z s p c w \ 

m z cpc w fi Cihnfi (2^22^ CaJiisA* 

dm z c w h\ x vcp 

( 2 m Z C w /l22 uc /8 ° 

V C3m z c w h e 33 vcp J 



(47) 



Note that all gaugino and higgsino mixings with neutrinos and charged leptons are 
parametrized by the three parameters Q. 

In the following section we shall need the couplings of gravitino, W- and Z-bosons 
to neutralino and chargino mass eigenstates. Since (j <1, diagonalization of the mass 
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Figure 1: The parameters \i and B of Eqs. |3J) and HJ), respectively, as functions of the bino mass Mi. 
The plot has been obtained by means of S0FTSUSY3 . J22J. 



matrices to first order in Q is obviously sufficient. We shall also consider supergravity 
models where the supersymmetry breaking parameters satisfy the inequalities (cf. Fig. [T]) 



m z < Mi )2 < fi ■ 



(4* 



The gaugino-higgsino mixings are 0(mz/fj), and therefore suppressed, and xl, the light- 
est neutralino, is bino-like. 

The mass matrices A4 N and Ai c are diagonalized by unitary and bi-unitary trans- 
formations, respectively, 



U^ T M N U^ = MZ ag , UWM C U® = Mg ag , 



(49) 



where £A n )t U^ = £/( c )t[/( c ) = [/( c )t£/( c ) = i. These unitary transformations relate the 
neutral and charged gauge eigenstates to the mass eigenstates (x°, v 1 ^) (a = 1, . . . , 4) and 
{Xai e 'di (Xa; e 'j C ) ( a — 1)2), respectively Inserting these transformations in Eqs. ( 1421) - 
(1441) and dropping prime superscripts, one obtains neutral, charged and supercurrents in 
the mass eigenstate basis: 



+ {tal^^i + X*l^f> e) e t + xll^^et + h.c.) , (50) 

: tai,v^ x ~ + xi^vjr^i + v a ,vt x) x- + vnMr^i > ( 5i ) 

S" =\ [Y, 1 P \ Y {u^ X °a + U^>x° a c + U^\ + U^ v Kt) F vp + . . . , (52) 



J, 
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where we have defined the photino matrix elements 

Up** = c w U^ + s w UJT^ , Up") = cjjf^ + s w ut' u) . (53) 

In the appendix the unitary transformations between gauge and mass eigenstates and the 
resulting matrix elements of neutral and charged currents are given to next-to-leading 
order in rriz/fi- As we shall see, that expansion converges remarkably well. 

In the next section we shall need the couplings of the lightest neutralino Xi to charged 
leptons and neutrinos, and the coupling of the gravitino to photon and neutrino. From 
the formulae in appendix A one easily obtainsj (S2/3 = IspCp) 

VS--0^(1 +0 (^)), (55) 

TT tf,v) r m z s w c w (M 2 - Mi) / / m|\\ 

Note that the charged and neutral current couplings agree up to the isospin factor at 
leading order in m|//i 2 , i.e., V^'^q = V^'l Q /2. The mass of the lightest neutralino is 
given by 

We have numerically checked that varying Mi between 120 and 500 GeV, the relative 
corrections in Eqs. (1511) - (157)1 are less than 10%. 



5 Fermi-LAT and the LHC 

We are now ready to evaluate the implications of recent Fermi-LAT data [T511T5] and 
cosmological constraints [23H25J for signatures of decaying dark matter at the LHC. We 
shall first discuss monochromatic gamma-rays produced by gravitino decays and then 
analyze the implications for a neutralino and a r-NLSP, respectively. 

In order to keep our analysis transparent we shall not study the most general pa- 
rameter space of softly broken supersymmetry, but only consider two typical boundary 
conditions for the supersymmetry breaking parameters of the MSSM at the grand unifi- 
cation scale, 

(A) m = mi/2, a = 0, tan/3 = 10 , (58) 



4 The matrix element U^' v agrees with the one used in [13,29] for M2 — Mi <C Mj.. 
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Figure 2: Gravitino decay into photon and neutrino. 

with equal universal scalar and gaugino masses, mo and m.1/2, respectively; in this case a 
bino-like neutralino is the NLSP. The second boundary condition corresponds to no-scale 
models or gaugino mediation, 

(B) m = 0, mi/2, ao = 0, tan/3 = 10 , (59) 

which yields the right-handed stau as NLSP. In both cases, the trilinear scalar coupling 
ao is put to zero for simplicity Choosing tan/3 = 10 as a representative value of the 
Higgs vacuum expectation values, only the gaugino mass parameter mi/2 remains as in- 
dependent variable; the mass parameters \x and B are determined by requiring radiative 
electroweak symmetry breaking with the chosen ratio tan/3. For both boundary con- 
ditions (158)) and (159~|) . the gaugino masses at the electroweak scale satisfy the familiar 
relations 

Mi M> 

-^ ~ 6.0 , -£- ~ 1.9 . (60) 

Mi Mi v ; 

For the chosen supergravity models, consistency with electroweak precision tests, 
gravitino dark matter (GDM) and thermal leptogenesis leads to the following allowed 
mass ranges of gravitino and lightest neutralino 



10 GeV < m 3/2 < 500 GeV , 100 GeV < m x <> < 500 GeV , (61) 

where we have used m x o ~ M x (cf. (|57[ )). Note that the masses Mi and m 3 / 2 cannot be 
chosen independently. The GDM constraint implies that for a given gravitino mass the 
maximal bino mass is JWf* ~ 270 GeV(m 3/2 /100 GeV) 1/2 [20]. 

Consider now the rate for gravitino decay into photon and neutrinqj [10J (cf. Fig. [2]), 

r 3/2W ._L£|^gz|, (6 2) 



^3/2{l 1/ ) denotes the sum of the decay rates into photon neutrino and photon antineutrino. 
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Inserting the matrix element (loT)]) one obtains the gravitino lifetime 



32^2 G F Ml M\Ml 



rrir 



^™=ifri&w=wy i+0 \'»i?))< (63) 

where a is the electromagnetic fine-structure constant, and we have defined 

C 2 = EC 2 - (64) 

i 

The corrections to the leading order expression in ([HSD are less than 10%. Using Eq. ([6"Uj) 
and M P = 2.4 x 10 18 GeV, one obtains 

Recent Fermi-LAT data yield for dark matter decaying into 2 photons the lower 
bound on the lifetime tdm(77) ^ 1 x 10 29 s, which holds for photon energies in the range 
30 GeV < E 1 < 200 GeV [18]. For gravitino decays into photon and neutrino this implies 

73/2 M >5x 10 28 s , 30 GeV < £ 7 < 200 GeV . (66) 

Since according to the GDM constraint the largest allowed bino mass scales like 
M™ ax ex m 3 , 2 , the largest lifetime ( |65|) . and therefore the most conservative bound on 
C, is obtained for the smallest value of m 3 / 2 - F° r small gravitino masses, a rough lower 
bound on the lifetime can be obtained from the isotropic diffuse gamma-ray flux. The 
recent Fermi-LAT data give E 2 dJ/dE\ 5 GeV ~3x 10~ 7 GeV (cm 2 s str) -1 [TS]. From the 
analysis in [12] one then obtains T3/2 > 10 28 sjj Together with Eq. ( 1651) one then obtains 
the approximate upper bound on the R-parity breaking parameter 

C < 3 x 10- 8 . (67) 

On the other hand, the observation of a photon line corresponding to a gravitino lifetime 
close to the present bound would determine the parameter £ aCl 

c - W 5 x l ° 28s \ 1/2 ( Mi U ^3 /2 v 3 / 2 rfi8 , 

Cobs " 10 \^MJ l200GevJ llbTGevJ (68) 

Note the strong dependence of Cobs on the gravitino mass. In (168D we have normalized 
these masses to central values suggested by thermal leptogenesis, electroweak precision 
tests and gravitino dark matter [2Uj . 
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Figure 3: Neutralino decays into charged lepton and W-boson, and neutrino and Z-boson. 

5.1 Neutralino NLSP 

A neutralino NLSP heavier than 100 GeV dominantly decays into charged lepton and 
W-boson or neutrino and Z-boson [30] (cf. Fig. [3]). The partial decay widths are given 
by 



G F m\ 



r ( x ; -, w*p) =1 -^ £ \V^& | Mm,;) (l + o ( S2 ^ 



(69) 
(70) 



Here V^ lq and V^lo are ^ ne charged and neutral current matrix elements at leading 
order, which are given in Eqs. (1551) and (1541) . respectively, and 



fw,z{m x o 



m 



w,z 



Xi> 



m 2 

XI 



1 + 2- 



m 



w,z 



rn? 

Xi 



(71) 



is a phase space factor which becomes important for neutralino masses close to the lower 
bound for m x? of 100 GeV (cf. Fig. H}. 

The total neutralino NLSP width is the sum 



r x? = r(x? -► w±n + r(xS -► ^) • 

Using the matrix elements (1541) and (1551) . one obtains the branching ratios 

BR (xl -»■ W^/ T ) ~ 2 Si? (xS -> ^) • 
Furthermore, the flavour structure of our model implies 

Si? (x? -> W^ T ) ~ BR (xi -> ^V 



(72) 



(73) 



(74) 



6 This lifetime is obtained by reseating in Fig. 2 of [12] the signal by the factor 0.1. 
7 The results ([57)1 and (|68[) are approximately consistent with the recent analysis [IT] . 



15 




100 150 200 250 300 350 400 450 



500 



m x o [GcV] 



Figure 4: Phase space suppression factor for neutralino decay to Z-boson and neutrino. 



Using the matrix elements (|54|) . (|55p and (|56|) for neutral, charged and supercurrent, 
respectively, one can express the neutralino lifetime directly in terms of the gravitino 
lifetime, 



T xl 



(M 2 - Mi) 2 m 



3/2 



Ml 



G F M 2 P m\ 

Al 



2\/2 

x (2/w(m x o) + fz(m x o 



^3/2 (7 



V) 



-1 




l + 0[ s 2B ^ 



(75) 



With the mass relations (I57p and (|6"0"j) one then obtains for the minimal neutralino decay 
length 



cr x o > 80 cm 



-3 



( m X° 1 \~ a ( ^3/2 \ 

U50GeV/ UOGeV/ 



x (2f w (m x o) + fz(m x o) 




1+ °'V, 



(76) 



In Eqs. (1751) and (1761) the corrections to the leading order expressions are less than 10%. 
We emphasize again the strong dependence of this lower bound on the neutralino and 
gravitino masses. For instance, for a gravitino mass of 100 GeV and the Fermi-LAT 
bound T3/2 > 5x 10 28 s, which applies for gravitino masses in the range 60 GeV < 
^3/2 < 400 GeV, one obtains ct x o ~ 2 km for a neutralino mass of 150 GeV. It is very 
interesting that such neutralino lifetimes are detectable at the LHC [31J. 

We conclude that, given the current bounds on the gravitino lifetime, a neutralino 
NLSP may still decay into gauge boson and lepton inside the detector, yielding a spec- 
tacular signature. However, for most of the parameter space a neutralino NLSP decays 
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outside the detector, leading to events indistinguishable from ordinary neutralino dark 
matter. 

5.2 r-Lepton NLSP 

Contrary to the neutralino NLSP decay, the R-parity violating decays of a ri-NLSP 
strongly depend on the flavour structure and the supersymmetry breaking parameters. 
The relative strength of the various decay modes becomes most transparent in the field 
basis where all bilinear R-parity breaking terms vanish, as discussed in Section 2. Since 
the R-parity breaking Yukawa couplings are proportional to the ordinary Yukawa cou- 
plings, decays into fermions of the second and third generation dominate. The leading 
partial decay widths of left- and right-handed r-leptons are (cf. (I14j) ) 

^r L (r R u) = ±-j: |A«3| V £ , (77) 



3 
lfor 1 



Tr L (hb R ) = T 7L (t L s R ) = TT-|A3 33 | 2 m~ L , (7* 



r ?L Ml) =1^1^3331 V*. (79) 

r~>^) = T~ R {ii L v) = ^ E IA.33I 2 m~ R . (80) 

In the flavour model discussed in Section 3, the order of magnitude of the various decay 
widths is determined by the power of the hierarchy parameter 77 (rf ~ 1/300), 

~ T~ L (t L b R ) ~ T~ L it L s R ) ~ V 4 Q 2 m 7 , (81) 

r~ L (t R b L ) ~ r) 2 e 2 m 7 . (82) 

The lightest mass eigenstate f± is a linear combination of tl and tr, 

T\ = sin r T£ + cos T f R . (83) 

From the above equations one obtains the Ti-decay width 

r~ = sin 2 6 T (T~ L (t r u) + 2r~ L (t L b R ) + T~ l (i R b L )) + 2 cos 2 6 t T 7r (t l u) . (84) 

The total width is dominated by the contributions f R — > tlP, \ii,v and tl — > i R bL, 
respectively, 

r~ = sin 2 6 T T~ L (t R b L ) + 2 cos 2 6 T Y~ R (t l u) , (85) 
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Figure 5: r -mixing angle: sin (8 T ) as function of the lightest r-mass m~ . 



and it can be directly expressed in terms of the r-lepton and top-quark masses, 



Tl 167TV 2 

where we have assumed 



(3m 2 t sin 2 6 T + 2m 2 . tan 2 /3 cos 2 6 T 



m 



n ' 



^2,3 



s 2,3 



: 2,3 



(86) 



(87) 



This corresponds to the parameter choice a = b = c=lin Eq. f)34p . Note that t\- 
decay width and branching ratios have a considerable uncertainty since these parameters 
depend on the unspecified mechanism of supersymmetry breaking. From Eqs. f)18p . (T2"6"j) 
and 7] ~ 0.06, one obtains for the R-parity breaking parameter 



e ~ C, — i]Q ~6x 10' 



(88) 



which is consistent with the present upper bound fIBTl) within the theoretical uncertain- 
ties. 

The dependence of the mixing angle 9 T on m~ is shown in Fig. for the boundary 
condition (I59p . For masses below the top-bottom threshold only leptonic rx-decays are 
possible. When the decay into top-bottom pairs becomes kinematically allowed, sin 2 9 T 
is small. However, the suppression by a small mixing angle is compensated by the larger 
Yukawa coupling compared to the leptonic decay mode. This is a direct consequence of 
the couplings A' which were not taken into account in previous analyses. 

Due to the competition between mixing angle suppression and hierarchical Yukawa 
couplings, the top-bottom threshold is clearly visible in the r^-decay length as well as 
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Figure 6: t\ -decay length as function of m~ . Above the top-bottom threshold hadronic decays decrease 
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the branching ratios into leptons and heavy quarks. This is illustrated in Figs. E] and 
El respectively, where these observables are plotted as functions of m~ . Representative 
values of the ri-decay lengths below and above the top-bottom threshold are 



m~ <m t + m h : ct~ x 



m~ > m t + m h : ct~ x 



= 1.4 m , (89) 

150 GeV V5 x 10- 8 / 

= 0.6mf — r) . (90) 

250 GeV V5 X lO" 8 ' 



Choosing for e the representative value (1551) from gravitino decay, e = Cobs = 10 -9 , one 
obtains cr~^ = 4 km(l km) for m~ = 150 GeV(250 GeV). It is remarkable that such 
lifetimes can be measured at the LHC 



Is it possible to avoid the severe constraint from gravitino decays on the fi-decay 
length? In principle, both observables are independent, and the unknown constants in the 
definition of e, e' and e" can be adjusted such that ( = 0. However, this corresponds to a 
strong fine-tuning, unrelated to an underlying symmetry. To illustrate this, consider the 
case where the soft R-parity breaking parameters vanish at the GUT scale, S, = m\ d = 0, 
which was discussed in Section 3. In bilinear R-parity breaking, also the R-parity violating 
Yukawa couplings vanish at the GUT scale. With the one-loop radiative corrections at 
the electroweak scale (cf. (13"9"j) ; e« = /ij//i), 

5,(Aew) = ^r~ 2 [\sf i M 1 + Qg 2 M 2 ) In ^1 m | d (A EW ) = , 
lo7r \5 / A EW 

and Mi ;2 ~ //, one reads off from Eqs. (TTU1) . (1T2|) and (1 13)) 

4< = ^fe)- (91) 

Hence, all R-parity breaking parameters are naturally of the same order, unless the fine- 
tuning also includes radiative corrections between the GUT scale and the electroweak 
scale. 

Even if one accepts the fine-tuning ( = 0, one still has to satisfy the cosmological 
bounds on R-parity violating couplings, which yield e^ = fii/fi < 10~ 6 [22] • In the flavour 
model discussed in Section 3 this corresponds to the choice a = 20 in Eq. (13"3"|) . For the 
smaller fi-mass, which is preferred by electroweak precision tests, one then obtains the 
lower bound on the decay length 



cr~ 



> 4 mm . (92) 

150 GeV V ' 



However, let us emphasize again that current constraints from Fermi-LAT on the diffuse 
gamma-ray spectrum indicate decay lengths several orders of magnitude larger. 
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5.3 Planck Mass Measurement 

It has been pointed out in |9J that, in principle, one can determine the Planck mass 
from decay properties of a r-NLSP together with the observation of a photon line in the 
diffuse gamma-ray flux, which is produced by gravitino decays. This is similar to the 
proposed microscopic determination of the Planck mass based on decays of very long 
lived f-NLSP's in the case of a stable gravitino |33j . 

From our analysis of NLSP decays in this section it is clear that neutralino NLSP 
decays are particularly well suited for a measurement of the Planck mass, which does 
not require any additional assumptions. Eq. (I75|) implies [Gp = \/2/(4f 2 )), 

M 2 \m x oJ \ t x o J 
x (2f w (m x o) + fz(m x o)y 1/2 (l + O U^)) • (93) 

As expected, for gravitino and neutralino masses of the same order of magnitude, the 
ratio of the two-body lifetimes is determined by the ratio of the electroweak scale and 
the Planck mass, 



T x1 v 

r ^ J 77 

r 3 / 2 (7z/) M P 



(94) 



Quantitatively, using the relation ( 1601) for the gaugino masses, one finally obtains (v 
174 GeV), 

M P =3.6 x 10- GeV ( ^^ ^^^ ' ™ " ^ 

\m x oJ 




9 

mi 



x (2f w (m x o) + f z (m x o)) f 1 + U^ J J ■ (95) 

It is remarkable that the observation of a photon line in the diffuse gamma-ray flux, 
together with a measurement of the neutralino lifetime at the LHC, can provide a mi- 
croscopic determination of the Planck mass. 



6 Summary and conclusions 

We have studied a supersymmetric extension of the Standard Model with small R-parity 
breaking related to spontaneous B — L breaking, which is consistent with primordial 
nucleosynthesis, thermal leptogenesis and gravitino dark matter. We have considered su- 
pergravity models with universal boundary conditions at the GUT scale, which lead to 
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scalar tau or bino-like neutralino as NLSP. Supersymmetry breaking terms have been 
introduced by means of higher- dimensional operators. The size of the soft terms corre- 
sponds to gravity or gaugino mediation. 

We have analyzed our model, which represents a special case of bilinear R-parity 
breaking, in a basis of scalar SU(2) doublets, where all bilinear terms vanish. In this basis 
one has R-parity violating Yukawa and gaugino couplings. They are given in terms of 
ordinary Yukawa couplings and 9 R-parity breaking parameters e,, e' t and e", i = 1, ..., 3, 
which are constrained by the flavour symmetry of the model. The R-parity violating 
couplings include terms proportional to the up-quark Yukawa couplings, which were not 
taken into accound in previous analyses. 

The main goal of this paper are the quantitative connection between gravitino decays 
and NLSP decays, and the corresponding implications of recent Fermi-LAT data on 
the isotropic diffuse gamma-ray flux for superparticle decays at the LHC To establish 
this connection one needs the relevant R-parity breaking matrix elements of neutral, 
charged and supercurrents. For the considered supergravity models these matrix elements 
can be obtained analytically to good approximation, since the diagonalization of the 
neutralino-neutrino and chargino-lepton mass matrices in powers of vrtzj [i converges 
well, as demonstrated in the appendix. The analytic expressions for the decay rates 
make the implications of the Fermi-LAT data for NLSP decays very transparent. 

Our main quantitative results are the branching ratios for NLSP decays and the 
lower bounds on their decay lengths. For a neutralino NLSP with m x o = 150 GeV, the 
Fermi-LAT data yield the lower bound cr x o > 30 cm. This bound does not depend 
on details of the superparticle mass spectrum or the flavour structure of the model. It 
directly follows from the comparison of two-particle gravitino and neutralino decays. On 
the contrary, there exists no model independent lower bound on the Ti-decay length. The 
natural relation between gravitino and r-decay widths can be avoided by fine-tuning. In 
this case the cosmological constraint that the baryon asymmetry is not washed out leads 
to the lower bound cr~^ > 4 mm. 

Without fine-tuning parameters the diffuse gamma-ray flux produced by gravitino 
decays constrains the lifetime of a neutralino as well as a r-NLSP. For typical masses, 
"^3/2 ~ 100 GeV and wnlsp ~ 150 GeV, the discovery of a photon line with an intensity 
close to the present Fermi-LAT limit would imply a decay length ctnlsp of several hun- 
dered meters. This is a definite prediction of a class of supergravity models. It is very 
interesting that such lifetimes can be measured at the LHC [3"T|I3"2]. 

Finally, it is intriguing that the observation of a photon line in the diffuse gamma- 
ray flux, together with a measurement of the neutralino lifetime at the LHC, can yield 
a microscopic determination of the Planck mass, a crucial test of local supersymmetry. 
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A Appendix: Gauge and mass eigenstates 

A.l Mass matrix diagonalization 

The mass matrices M. N and Ai c in the gauge eigenbasis were explicitly given in Eqs. 
and (jUD , respectively, 



M 



N 



( Mi 


mzSpSu, 



m z sps w -m z C/3S w 

M 2 -m z spc w m z cpc w 

-m z spc w —fi 



M 



c 



-m z cps w m z cpc w -fi 

V -dmzs w Qm z c w 

/ M 2 m z s p c w \ 

m z c p c w n Ci^iiA* (2h e 22 n (3^33^ 

dm z c w h\ x vc$ 

(2mzc w h% 2 vcp 

V C,?,m z c w h% z vcp J 

For non- vanishing i?-parity breaking parameters Q, i = 1, 
between gauginos, Higgsinos and leptons, 



-dfn z s w \ 
dm z c w 





i 



3, they induce a mixing 



-C D^ (b, w 3 , hi h° d , V% ) M N (b, w\ hi hi v % 



+ ({w~,K, ei )M c (w + } hlef) 



+ h.c. 



(A.l) 



The matrices A4 N and M. c are diagonalized by unitary and bi-unitary transformations, 
respectively, 



U (n)T M N u (n) = M g^ ^ V ^ M ° U^ = Mg ag , 



(A.2) 



where [/Mt £/("■) = £/( c )t[/( c ) — £/( c )t[/(c) _ i These unitary transformations relate the 
neutral and charged gauge eigenstates to the mass eigenstates (xl v 'i) (a = 1, . . . , 4) and 
(Xai e 'i)i (Xa,ef) (a = 1,2), respectively. 
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In this work we consider the two boundary conditions (A) and (B), defined in 
Eqs. (158)) and (1551) . respectively. The corresponding supergravity models satisfy the re- 
lation (HED, m z < Mi )2 < /i, and in the regime 120 GeV < M 1 < 500 GeV one finds 
0.07 < vnzl ' \x < 0.25. We diagonalized the above mass matrices to first order in the small 
parameters Q and to second order in mz/ [i- The size of the relative corrections given 
below has been calculated for the above parameter range using S0FTSUSY3.0 [23]. As we 
shall see, the relative corrections are of order m|//z 2 , and the expansion converges well 
for most matrix elements. 



The neutralino and neutrino mass eigenvalues are 






A2 



X3 



m v o 



m h 



M l - 



M 2 



/' 



/i 2 - Ml 
m\ (M 2 + fjLs 2 p) 



c„ 



V? - M 2 




m z (fi-M 1 c 2 w -M 2 sl)(l + s 2f} ) 
2 (ji - MO (// - M 2 ) 

m|(/i + M 1 4 + M 2 4)(l- S2/3 ) 




/i 



o + o k 



2 (// + Mi) (// + M 2 ) 



.'^Kf, 



2 m | N 



(A.3) 
(A.4) 
(A.5) 
(A.6) 

(A.7) 



We checked numerically that relative corrections 0(m 2 z /fi 2 ) to the above neutralino 
masses are smaller than 0.05, 0.15, 0.10, 0.001, for m x o, . . . , m x o, respectively. 

The chargino and lepton mass eigenvalues are 



m xt 



m xf 



m P > 



Mo 



/i + 



m z 2 (M 2 + fis 2f3 )c. 



»i 



tf 



Mo 



m z 2 {\x + M 2 s 2 p)c. 



2(/i 2 



Mo 2 



2 m Z 



h&cp I 1 + o I c //2 




(A8) 

(A.9) 

(A.10) 



Here the relative corrections of 0(m 2 z /fi 2 ) are numerically smaller than 
The unitary matrix [/( n ) from Eq. (1A.2I) can be written as 



/o. 



lj(n) 



U, 



(x°) 
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ab 



U, 



(xV) 



u, 
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(A.ll) 
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with 
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+ 



^ 1 o \ 

10 





m z 2 (Mi 2 +2/iS 2/3 Mi+^ 2 )s 



1 1 

— — 

72 y/2 I 



2{M 1 2 - f i 2 ) 

rn z 2 (Mi+iJ,S20)s2w 

2{M 1 -M 2 )(M 1 2 - IJ ?) 

mz(lJ,cp+Misp)s w 

Mi 2 -^ 2 
rn z (M 1 cp+^s /3 )s w 



m z 2 {M2+lJ,S2f i )s2w 
2(Mi-M 2 )(M 2 2 -/i 2 ) 

m z 2 c 2 ,(M 2 2 +2^ 2 gAf 2 +M 2 ) 

2(M 2 2 -m 2 )" 
m z c w ([j,cp+M2Sp) 



m z (c l3 +s l3 )s l 



m z {c l3 -s l3 )s w \ 
V2{Mi+v) 



m z c w (c0+sp) m z c w (sp-c l3 ) 
V2(M 2 -ij) V2(M 2 +fi) 



M 2 2 - M 2 



M{ 2 -fi 2 



m z c w (M2Cp+iisp) 
M 2 2 - M 2 



m z 2 (cp+sp) 



Xi 

J 

-x 3 



{c0-s )m z 2 

J? X 2 

(cp-sp)mz 
/,2 



o 



rrir 



F 



where we used the abbreviations 

jj_ / (M 2 g /3 -(M 2 -2/i)c /3 )4 (M lg/3 - (Mx - 2/i) C/3 )4/ 
4^1 (M 2 -/i) 2 + (M^/i) 2 
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x 4 y 
(A.12) 

(A.13) 
(A.14) 
(A.15) 
(A.16) 



(M 2 + /x) 2 (Mx+^u) 2 

The numerical error of the matrix (IA.12I) in our parameter range of interest is smaller 
than 40% of the given NLO term. We do not discuss the slow convergence for this Im- 
parity conserving sub-matrix further, since this is beyond the scope of our analysis. 

Furthermore, 
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(A.17) 
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(A.18) 



(A.19) 
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The uncertainties in Eq. ( 1A.17I) evaluate numerically to less than 5%. For U^' x they 
are less than 0.15, 0.10, 0.25, 0.25, for a = 1, . . . , 4, respectively. 

The unitary matrices U^ and U^ c ' which diagonalize the matrix Ai c , cf. Eq. (1A.2I) . 
can be denoted as 
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(A.20) 



We find 
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(A.21) 

(A.22) 



Numerically, the relative correction to the NLO contribution to U^i is less than 25%. 
The off-diagonal elements of the matrix U^ c ' to leading order in h^ are 



U, 



(x + ,e c ) 







1 



m z c w (M2Sfj-vcp) 

J- t h e I *^ 

r ^ n ii \ m z 2 c l3 cl,(v l ic l3 +M2(v-Li)sf3) 

M 2 V 



u, 



(e c ,X + ) 



CM 



o 



32/3"^ 
M 2 



(A.23) 



+ CA e 



M 2 2 (M|- M 2) 
m z 2 /j 2 cg,^ 



2(y-M 2 2 )' 



1 + 



ml 



/.;/ 



(A.24) 



where 



1/ =— (^s 2/3 M 2 3 + /i (2t;M 2 2 + /i 2 (/i - 2v)) c 2 + lisp U (2/i 2 - M 2 2 ) S/9 
/i 



-2M 2 (M 2 2 + (v- 2//)// 



c /3 



(A.25) 



The numerical relative correction to the NLO term in U^ a ' x is smaller than 0.10,0.15 
for a = 1,2, respectively. For U[f ' e it is smallecl than 1%, and smaller than 10% for 

fr(x + ,e c ) 
<^2i 



8 The numerical calculation of the error reaches our numerical precision. The given value is calculated 
from the comparison with the analytical NNLO expression. 
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The block diagonal elements of the matrix U^ are 



U, 



Or) 



a/3 



1 
1 



/ m z 2 c^(M 2 c / 3+/iS /3 ) 2 m z c w (M2Cp+nsp) \ 



2(M 2 2 - M 2 ) 

m z Ct 1 ,(Af 2 c^+M s /3) 
Af 2 2 -^ 2 



mz 2 cS.(^2C^+/j^) 2 
2(M 2 2 - M 2 ) 2 J 



l + O 



rrir 



P 



Ulf = 8 ij + Q (C 2 ) 



(A.26) 

(A.27) 



Numerically, the relative correction to the NLO contribution to U^i is smaller than 20%. 
The off-diagonal elements of U^ are 



u (x ,ej = c- 



m z c m 

M 2 
m z 2 cj 1 sp 

M 2 n 



l + O 



u ia Si 



m z c w 

M 2 

m z 2 c 2 u (ticp+M2Sp) 

ti 3 -M 2 2 ii 



g2/3^| 



l + O 






(A.28) 
(A.29) 



Here we ignored corrections that are proportional to the Yukawa couplings h\ { or higher 
powers thereof. The numerical value of the higher order correction relative to the NLO 
term is smaller than 1% for XJ X , smaller than 5% for U^ ,x , and smaller than 15% 
for U$ x ~\ 

A. 2 The currents in mass eigenstate basis 

The neutral and charged currents were given in Eqs. floU]) and (J5T|) . 



-— 0„, T/(x°),,0 i — - T /(X L- i ^t+„, T/(x + )„+ i „„, T/(") T , i --. T^( e ). 



^ =x»7mK6 j x" + x*i»v$ % + X*1»V$ \% + VaXi ^ + ^7/,Ki '< .- 



(X .e)„ _+ T/ (x + ,e<+c 



+ (x^K?'^ + 5£7|J# ~' 6) e* + X+7mK 



e- + h.c. 



J, = tal^ha+tal^ X ' e) e i +Va^Vt X \-+Va^' e) e j . 

The CKM-like matrices vg\ V^\ V$ + \ V$\ if, V^\ V^f' e \ v£ + > e °\ V<£, 
Val ■> Via' i Vij follow from the currents in gauge eigenbasis, Eqs. (jl2"j) and (J13j) . and 



the explicit matrices U^ n \ U^ c > and U^ c \ Here we focus on the matrices V^; and V^, 
since they determine the interactions of interest for this work. We find 

/ s w m z \ ( I s 20 m 2 z \^ 

l + O 



v v 



(X,v) 







s w m z 

2A'h 
c w m z 
2M 2 
m z 2 fivi 



2v / 2Mi(Mi-/j)(A'/2-m) 

m z 2 nv 2 

\ 2V2M 1 (M 1 +f J ,)(M 2 +n) ) \ 



\ 



m z 
m z 

»\ 

?2pmt 



' ai 



(A.30) 
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with abbreviations 



Vl 



V-2 



Mi (Mi - fj)((fi - 2M 2 ) C/3 - fisp)c: 



M 2 /i 2 

+M 2 (M -/')((/' -2M l )r-,.-/,.s ; ).s 
1 



2 

ui y ? 



Mi (Mi + /i)((2M 2 + /i) C/3 + /i S/3 K 



M 2 /x 2 

+M 2 (M 2 + a*)((2Mi + /i) C/3 + ^4 



/ 



V, 



(X,e) 



S m m z 



\ 



Mi 
(v / 2-l)c tu m, z 

M^ 

M 2 (Mi-n}(M 2 -fJ.) 
V M 2 (Mi+ m )(M 2 +m) / 



/ 



/ s 2l3 m 2 z \ \ 



1 + 



v i# y 



with abbreviations 

1 — 

«i =2~2 ( c /3 (( M i " M) (v^M, - 2/i) 4 + v / 2M 2 (M 2 - /j,)s\ 

+ S/3 ((72 - 2) (Mi - fi)fic 2 w + v / 2M 2 (M 2 - fi)s 2 w )) , 

V2 =-^ (cfi ((Mi + /i) (v^M, + 2») 4 + v / 2M 2 (M 2 + ^ 



2/x 2 



-S/3 



V^ - 2) fi(M 1 + /i)4 - v / 2M 2 (M 2 + ,u)< 



(A.31) 



(A.32) 



Numerically, the relative errors are smaller than 0.10, 0.20, 0.15, 0.05 for a = 1, . . . , 4. 
Finally, 



(A.33) 



(A.34) 



(A.35) 



Here we again neglected corrections that involve the Yukawa couplings h e i{ . The numerical 
corrections to the NLO contributions to V^f are smaller than 0.05,0.15,0.20 for a = 
1, 2, 3, respectively. For a = 4 we reach the limit of our numerical precision. 
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